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ABSTRACT 

We present Submillimeter Array observations toward the 10^'^ L0 star- 
forming region G240.3H-0.07, in the J = 2-1 transition of ^^CO and ^^CO and at 

1.3 mm continuum, as well as the ^^CO and ^^CO observations from the Caltech 
Submillimeter Observatory to recover the extended emission filtered out by the 
interferometer. Maps of the ^^CO and ^^CO emission show a bipolar, wide-angle, 
quasi-parabolic molecular outflow, roughly coincident with an IR nebula revealed 
by the Spitzer 3.6 and 4.5 fim emission. The outflow has ~ 98 Mq molecular 
gas, making it one of the most massive molecular outflows known, and resulting 
in a very high mass-loss rate of 4.1 x 10^^ MqJt~^ over a dynamical timescale of 

2.4 X 10^ yr. The 1.3 mm continuum observations with a 4" x 3" beam reveal a 
flattened dusty envelope of ~ 150 Mq, which is further resolved with a 1.2" x 1" 
beam into three dense cores with a total mass of ~ 40 Mq. The central mm 
core, showing evidence of active star formation, approximately coincides with 
the geometric center of the bipolar outflow thus most likely harbors the powering 
source of the outflow. Overall our observations provide the best case to date of 
a well-defined wide-angle molecular outflow in a > 10^ star-forming region. 
The outflow is morphologically and kinematically similar to low-mass protostellar 
outflows but has two to three orders of magnitude greater mass, momentum, and 
energy, and is apparently driven by an underlying wide-angle wind, hence further 
supports that high-mass stars up to late-0 types, even in a crowded clustering 
environment, can form as a scaled-up version of low-mass star formation. 

Subject headings: ISM: individual (G240.31+0.07) — ISM: jets and outflows — 
stars: formation — stars: early-type 
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Introduction 



Single-dish surv eys have shown that molecular outflows are ubiquitous in high-mas s 
star forming regions (IShepherd &: Churchwelllll996l : IZhang et al.ll200ll : iBeuther et al.ll2002al ). 
According to the standard paradigm of low-mass star formation this provides a statistical 
argument that high-mass stars may form via accretion processes. However, given the typical 
large distances (> 1 kpc) and clustering environment of massive star formation, single-dish 
observations of > 10" resolutions cannot resolve the detailed structures of massive molecular 
outflows or identify their driving sources. Recently there have been an increasing number 
of interferometric studies which p rovide new insig hts into the morphology and kinematics of 
massive molecular outflows (e.g.. iQiu et al.ll2007l . and references therein), but the statistics 
is still poor and many basic questions, e.g., how massive molecular outflows systematically 
differ from their low-mass counterparts, how they are driven, and how they affect their 
surroundings and host clouds, remain open. To address these questions, which are essential 
to understand massive star formation, a large set of high-angular-resolution observations 
covering a wide range of luminosity and evolutionary stages is needed. 

Toward these goals we study G240.31+0.07 (hereafter G240), an active hi gh-mass star- 
form ing region with a far-infrared luminosity of 10^'^ Lq at a distance of 6.4 kpc (IMacLeod et al 



19981). 



t harbors an ultra-compact (UC) Hii region and is associated with H^O and OH 



Kumar et al 



masers (IHughes fc MacLeodlll993l : ICaswell Ill997l : MacLeod et al.lll998l : Migenes et al.lll999r ) . 



J2003h reported a parsec-scale velocity gradien t in C^^O (2-1) wit h a 33" beam 



and the detection of near-infrared H2 emission in this region. I Chen et al.l (120071 ) presented an 
interferometric study at 654 GHz of the central core and resolved the dust continuum emission 



pes: 


SheD 


lerd & Churchwell 


iHunter 


( 


1997 


) mapped the 



CO (3-2) emission with a 20" beam and proposed two overlapping outflows both in the 
northwest-southeast orientation. However, limited by low angular resolutions, neither of the 
two single-dish studies provided a detailed analysis of the outflow morphology and kine- 
matics. High-angular-resolution observations are crucial to unveil the outflow nature of this 
luminous source. 

Here we present a ^^CO, ^^CO (2-1) study of the G240 outflow, using the Submillimeter 
ArraJ (SMA) observation. Complementary short spacings are obtained from the Caltech 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
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Sub millimeter Observator}o (CSO) 10.4 m telescope. The archival Spitzer mid-infrared 
imaging is also presented for comparison. The 1.3 mm continuum observations are used to 
pinpoint the outflow powering source. 



2. Observations 

The SMA observations were carried out on 2008 February 25 with 8 antennas in the 
compact configuration and on 2008 February 16 with 7 antennas in the extended configu- 
ration. The 2 X 2GHz correlator was configured to cover ^^CO (2-1) in the upper sideband 
and ^^CO (2-1) in the lower sideband. And the continuum was constructed from line-free 
channels. Lines of ^^CO, ^^CO (2-1) and the 1.3 mm continuum, which are relevant to 
the outflow study, are presented in this work. Many other lines covered in the same setup 
will be presented in a future paper. To ensure the coverage of the entire outflow we ob- 
served two flelds centered at (R.A., decl.)j2ooo = (07M4'°52.49^ -24^07"^52.1^) and (R.A., 
decl.)j2ooo = (07M4"^51.07^ -24^07""34.9"), respectively. We used Titan as the primary flux 
calibrator and 3c273 as the bandpass calibrator. The time dependent gain was monitored by 
observing quasars J0730-116 and J0826-225 every 23 mins. Visibilities were calibrated using 
the IDL MIR package and then output to MIRIAD for imaging. With natural weighting 
the synthesized beams in the compact and extended conflgurations are about 4" x 3" and 
1.2" X 1", respectively. 

The shortest baseline in our SMA observations is 16.5 m, corresponding to a spatial 
scale of 20" for observations at 225 GHz. Thus spatial structures more extended than 20" 
were not sampled in the SMA observations. This spatial flltering can signiflcantly affect 
the ^^CO, ^^CO (2-1) maps at velocities close to the cloud velocity. To recover the missing 
short spacing information we observed the ^^CO, ^'^CO (2-1) emission with the CSO 10.4 m 
telescope on 2008 February 12. During the observation the weather condition was excellent 
for 1 mm waveband with 'r225GHz ~ 0.08. The observations were carried out in the on-the-fly 
mode centered on (R.A., decl.)j2ooo = (07^^44^52. l^ -24'i07°>49"). We obtained a 15 x 15 
grid in ^^CO and a 13 x 13 grid in ^'^CO, with an integration of about 5 s on each cell. 
The 10" cell spacing is about | of the CSO beam, which is ~ 32.5" in ^^CO (2-1). By 
observing Mars and Saturn we derived a beam efficiency of 0.58 ± 0.10. The spectrometer 
used has 1024 channels throughout the 50 MHz bandwidth, resulting in a spectral resolution 
of 0.0488 MHz (about 0.065 km s~^) per channel. The final maps were smoothed to 2 km 
s~^ per channel. The data were reduced using the standard CLASS package. We combined 
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the SMA compact and CSO ^^CO, 
Zhang et all JlQQsh . 



13 



CO data in MIRIAD following a procedure outlined in 



The deep Spitzer IRAC images were obtained from the Spitzer archive (PID: 20499). We 
found no obvious artifacts from a visual inspection, thus adopted sl4.0 Post Basic Calibrated 
Data products provided by the Spitzer Science Center. 



3. Results 

3.1. mm Continuum Emission 

In Figure [T^, the 1.3 mm continuum emission observed with the SMA compact array 
shows an elongated envelope with its major axis in the northeast-southwest direction, as well 
as a minor peak to the northeast. At a higher angular resolution of 1.2" x 1" the central 
part of the structure is resolved into three dusty cores, namely MMl, MM2, and MM3 from 
north to south (Figure [T)d). Assuming optically thin dust emission, the dust mass can be 
estimated according to 

Mdust 



where M^ust is the dust mass, the continuum flux at frequency D the source dis- 
tance, ByiTdust) the Plan ck function at dus t temperature T^usti and = 10(z//1.2 THz)^ in 



cm^g~^ the dust o pacity (IHildebrand 1119831 ). Based on SMA observations toward this region 
Chen et al.l (120071 ) derived a gas rotational temperature of 96 K by comparing the brightness 
of two H2CO transitions, and derived an opacity index /? of 1.5 by comparing flux densities 
at 223 and 654 GHz. We then adopt Tdust = 100 K, (3 = 1.5, and a gas-to-dust ratio of 100, 
to estimate the mass of the larger envelope and the combined mass of the three dusty cores, 
which amount to 150 Mq and 40 Mq, respectively. The uncertainty of this estimate mainly 
comes from the determinations of (3 and T^usu for example, deviations of ±0.3 in j3 or ±50 
K in Tdust will induce ~ 50% uncertainty in the mass estimates. 



3.2. Bipolar Molecular Outflow 

3.2.1. Outflow morphology and kinematics 

Figure [2] shows the integrated ^^CO, ^^CO (2-1) emission observed with the SMA, as 
well as the Spitzer IRAC imaging for comparison. Because the interferometric filtering to 
extended emission helps to distinguish the outflow from ambient gas, we present the SMA 
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observations without merging with the single-dish data. The cloud velocity [vdoud) with 
respect to the local standar d of rest is about 67 .5 km s~^, which is adopted from the CSO 



C^^O (2-1) observations by iKumar et all ((20031) and consistent with our SMA C^^O (2-1) 



observations. The low- velocity channels, i.e., 58-64 km s~^ for the blueshifted lobe and 74-80 
km s~^ for the redshifted lobe in Figure [2^, are chosen to highlight wide-angle structures and 
to avoid contamination from diffuse gas in the surrounding, and the high-velocity channels 
(Figure [2]d) are those in the outer line wings. The low- velocity ^^CO emission clearly reveals 
a northwest-southeast (NW-SE) bipolar outflow centered at the group of the mm peaks 
and extending about 1.3 pc from end to end (Figure [2^). The outflow shell highlighted 
by limb-brightening outlines a quasi-parabolic structure with a wide opening angle. Both 
the spatial separation of the two lobes of such a wide-angle outflow and the partial overlap 
of blue- and redshifted emission in the NW lobe suggest that the flow axis is moderately 
inclined to the plane of sky. We suggest the molecular outflow seen in the ^^CO emission 
as ambient gas swept-up by an underlying wide-angle wind (see §4.20 . In this context, the 
NW wind is overall moving away from us, sweeping up ambient gas to create the redshifted 
^^CO outflow, but could be partly moving toward us in its near side, causing the blueshifted 
^^CO emission. The patchy appearance of the blueshifted emission might be attributed 
to the localized density structure of the medium. The SE counterpart of the wind may be 
expanding into a medium which is less dense in the far-side and would not create appreciable 
redshifted ^^CO emission. Immediately to the southeast of the bipolar outflow, there is a 
minor extension with overlapping red- and blueshifted emission, which might arise from 
another outflow. Meanwhile, the high-velocity ^^CO emission coincides with the inner part 
of the low-velocity bipolar outflow, with the bulk blueshifted emission revealing a clump 
in the middle of the low-velocity SE lobe and the redshifted emission a conical structure 
narrower than the low- velocity NW lobe (Figure [2}d). The bipolar outflow is also detected 
in the ^^CO emission with a morphology similar to but slightly less extended than the ^^CO 
outflow (Figure (St). 

Spitzer IRAC images often reveal outflow cavities through scatted light in the 3.6 /xm 
band and molecular outflow drivin g agents (e . g., un derlying jets or winds) through shocked 



H2 emission in the 4.5 /im band (iQiu et al.l l2008l ). In Figure the IRAC three-color 



composite image of the G240 central area yields a NW-SE nebula in the 3.6 and 4.5 fim 
emission. The north to south shell of the nebula appears prominent in the 3.6 fim band 
(blue) and coincides well with the north to south edge of the ^^CO outflow, suggesting its 
emission being dominated by scattered light from the outflow cavity wall, whereas the nebula 
structure coincident with the main NW-SE ^^CO lobe shows diffuse 3.6 and 4.5 /xm emission, 
which might be dominated by scattered light and shocked H2 emission. 

Figures [3] and m show channel maps of the ^^CO emission from the SMA observations 
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alone and from the combined SMA and CSO data at low velocities, respectively. It is 
immediately evident that for channels around Vdoud the combined data effectively recover 
the extended emission filtered out by the interferometer and significantly improve the side- 
lobes of the interferometric observations. In particular the diffuse gas emissions at 66 and 68 
km s~^ are mostly filtered out by the interferometer and recovered with the combined data. 
From a comparison among the spectra constructed from the CSO, SMA, and the combined 
data (Figure [5]), the SMA compact array observations recover less than 10% of the CSO flux 
at the line center, and about 80% of the CSO flux at the line wing of 54 and 80 km s~^, 
whereas the combined SMA and CSO data recover more than 85% of the CSO flux all the 
way from 54 to 80 km s~^; at the outer line wing of <52 km s~^ and >82 km s~^, there is 
little missing flux in the SMA observations. 

Outflow features emerge at 64 and 70 km s~^; the SE wide-angle shell at 64 km s~^ 
appears less prominent in Figure H] than that in Figure [3] due to the confusion from diffuse 
gas. Emissions from outflow structures appear predominant in channels of <62 km s~^ and 
>72 km s~^. The two lobes of the wide-angle outflow accord with each other in orientation 
and have similar opening angles at the base. They both appear more collimated at higher 
velocities: the NW lobe has a fan-shaped morphology at low velocities (most prominent at 72 
and 74 km s~^) and a narrower, F-shaped appearance at high velocities (most prominent at 
84 to 88 km s~^); the SE lobe shows the wide-angle shell at low velocities (most prominent at 
62 and 64 km s^^) and a clump away from the mm peaks at high velocities (most prominent at 
48 to 52 km s^^). Again the difference in detailed velocity structures of the two lobes could be 
at least in part attributed to the localized density structure of the environment within which 
the underlying wind is expanding. As seen in the integrated emission, a possible secondary 
outflow to the southeast of the mm peaks can be found in both blue- and redshifted channels. 

Channel maps of the ^^CO emission are shown in Figures Eland [71 for the SMA alone and 
the combined SMA and CSO data around Vdoud, respectively. Channels of 64-70 km s~^ are 
dominated by surrounding diffuse gas but also show extension along the outflow orientation. 
There is an embedded dense clump showing a northeast-southwest velocity gradient, which 
is also detected in other dense gas tracers. Emission from the bipolar outflow dominates 
channel of <62 km s~^ and >72 km s~^. The ^^CO emission appears confined within 10 km 
s~^ from Vcioud, for both blue- and redshifted lobes. 
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3.2.2. Outflow mass and energetics 



We use the combined SMA and CSO data to estimate the gas mass of the outflow, 
following 

Ti2 



M. 



16 59 

out - 1.39 X 10-'^exp(— ^)(Te, + 0.92)0^ 



where Mout, Tex, D, T12, and Si, are the outflow gas mass in Mq, excitation temperature of 
the ^^CO (2-1) transition, source distance in kpc, optical depth for ^^CO (2-1), and line flux 
in Jy, respectively. We adopt an excitation temperature of 30 K, a ^^CO-to-H2 abundance 
of 10~^, and a mean gas atomic weight of 1.36. Assuming the same excitation temperature 
for ^^CO (2-1) and ^^CO (2-1), we calculate according to 

^^(^^CO) _ 1 - e-"i^ 
^^(i3C0) ~ 1 - ' 

where x is the abundance ratio of ^^CO to ^"^CO, which is est imated to be 103 at a Galac- 
tocentric distance of 12.7 kpc for G240 (jwilson fc Roodlll994l ). Many ^^CO outflow studies 
estimate outflow mass either by assuming that ^^CO emission is optically thin at all outflow 
velocities or by correcting for a constant ^^CO opacity. Recently a few interferometric out- 
flow studie s use single-dish ""^^CO/^^CO line ratios to correct for ^^CO opacity as a function 
of velocity (jShepherd et al.lll998l . l2003l : ISu et al.ll2004l ). Our interferometric, beam-matching 
observations of ^^CO and ^^CO emission, both complemented by short spacing information 
from single-dish observations, enable a measurement of spatial and spectral distribution of 
^^CO/^^CO line ratios, thus a more sophisticated correction for the ^^CO opacity. We as- 
sume the ^^CO emission is optically thin where the ^^CO is not detected. The momentum 
and energy of the outflow is then derived following 



Pr 



out 



out 



V]V 



and 



E, 



out 



-^Mout{v)v^ 



respectively, where (and hereafter) v denotes the outflow velocity relative to Vdoud- To 
avoid contamination from the ambient molecular gas, we use channels of < 62 km s~^ in the 
blueshifted outflow and channels of > 72 km s~^ in the redshifted outflow, although channels 
of 64 and 70 km s~^ show contributions from the outflow (see Figures [3] and H]) . The area 
taken into account is defined by a polygon encompassing the outflow structure shown in 
Figure [2^. The dynamical timescale of the outflow is derived following tdyn = Lfio^j/vmax, 
where Ljiow ~ 0.65 pc and Vmax ~ 26.5 km s~^ are the half length of the end-to-end flow 
extension and the maximum flow velocity relative to Vdoud, respectively. Consequently the 
mass loss rate of the outflow can be derived following M^ut = Mout/tdyn- The derived outflow 
parameters are listed in Table [H 
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Discussion 



4.1. The outflow powering source 



From Figure [21 it is evident that one of the three mm continuum cores harbors the 
driving so urce of the bipolar outflow . From the maser and cm continuum observations in th e 
literature Jflughes fc MacLeodlEoQsl : ICasweU IEqQtI : iMigenes et aPllQQol : iKumar et allboosh , 
MMl roughly coincides with an UC Hii region and a H2O maser; MM2 is associated with H2O 
and OH masers but shows no cm emission down to 0.5 mJy at 6 cm or 2 mJy at 3.6 cm; MM3 
has no apparent signature of active star formation. The 1.3 mm continuum emission from 
MM3 appears much weaker than that from MMl and MM2, implying a mass signiflcantly 
less than 10 Mq. Its relatively faint dust emission and the apparent non-association with 
active star formation suggest that MM3 does not harbor the driving source of the bipolar 
outflow. Of the remaining two mm cores, MM2 appears closer to the geometric center of the 
bipolar outflow, thus is more likely to be the parent core of the outflow driving source. This is 
compatible with the relatively short outflow dynamical timescale, 2.4 x 10'^ yr, in comparison 



with > 10^ yr for massive molecular outflows from UC Hii regions (IShepherd et al 
2003h . 



1998 



4.2. Comparison with low-mass molecular outflows 



Low-mass molecular outflows are known to exhibit highly coUimated jetlike structures 
and/or wide-angle quasi-conical or quasi-parabolic structures. It is still unclear whether 
massive molecular outflows have an essentially similar appearance as low-mass outflows. 
High-angular-resolution interferometric observations have resolved a few single-dish bipolar 
outflows into multiple jetlike outflows, which were suggest ed as indirect eviden c e that as- 



socia t ed high-mass stars may form via disk accretion (e.g., iBeuther et al.ll2002bl : lOiu et al. 



2OO7I : Izhang et all l2007al jbl) . On the other hand, very few well- deflned wide-angl e molecu- 
lar outflows have been found in high-mass star-forming regions. lOiu et al.l (120081 ) reported 
bipolar conical outflow cavities in Spitzer mid-IR imaging of IRAS 20126+4104 and HH 
80/81, but the former is subject to ambiguity about whether it is caused by a precessing 
jet and the latter lacks kinematical information due to the lack of high-angular-resolution 
mm line observations. The G192.16-3.82 outflow, originating from an UC Hii region, may 
be representative of a bipolar wide-angle outflow with limb-brightened shells in a 10^ Lq 
star forming region (IShepherd et al.l 119981 ). Our observations of the G240 outflow reveal a 
bipolar wide-angle molecular outflow with a quasi-parabolic shell. To our knowledge this 
is hitherto the best case of a well-deflned wide-angle molecular outflow in a > 10*^ star 
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forming region. The outflow fia s a morpfiology simila r to protostellar low-mass niolecu- 



lar outflows, e.g L12 28 outflow flArce fc Sargentll2004f ). B335 outflow flHirano et al.lll988 



J0rgensen et al.l 120071 ). The behavior of higher collimation at higher velocities is also often 
seen in low-mass molecular outflows. 

Given the striking similarity between the G240 outflow and low-mass wide-angle out- 
flows, the G240 outflow could be driven by an underlying wide -angle wind, as commonly 



suggested for low-mass wide-angle outflows. IChurchwell I (119971 ) argued that neither accu 



mulated stellar winds nor entrained interstellar matter in bipolar jets is able to account for 
> 10 Mq of gas masses measured in massive molecular outflows typically at ~ 10^ yrs old. It 
is true that the G240 molecular outflow is unlikely to primarily contain accumulated stellar 
winds, because stellar winds from OB stars are radiation driven and their momentum supply 
rates cannot be higher than Lhoi/c, while the G240 outflow has a mechanical force, Pout/tdyn, 



by 


Churchwell 


1991 


). In this 



mixing layer. Because the environment cannot react to the rarefaction by the entrainment 
at velocities faster than the sound speed Cq (typically an order of 1 km s~^), the entrainment 
rate per unit area is given by epoCo, where e <^ 1 is the entrainme nt efficiency and po the 



density of the ambient medium. On the other hand, IShu et al.l (Il99ll ) proposed an analytical 
model for bipolar molecular outflows, where an underlying wide-angle wind radially sweeps 
up ambient gas into an expanding shell to form a molecular outflow. Further analytical 
calculations and simulations of the model, by adopting an angle-dependent wind blowing 
i nto an angle-dependent molecular cloud core, app ear generally consistent with observations 



(Li & Shu 


1996; 


Lee et al. 


2001; 


Shang et al. 


2006) 



sumed to be momentum driven and the entrainment rate per unit area is simply proportional 
to poVs, where Vg is the shell expanding speed (i.e., the molecular outflow speed observed in 
^^CO, typically an order of 10 km s~^) and determined by the momentum supply rate of the 
underlying wind. The mass loading limit in the turbulent entrainment model does not seem 
to be applicable to the swept-up shell model. Even in the turbulent entrainment model, 
a factor of 10 higher mass loading could be achieved if the entrainment process works at 
Alfven speeds, which are a factor of 10 higher than cq, when the magnetic flelds are impor- 
tant. Another molecular outfl ow model that can readily account f or high mass loading is the 



( Churchwell 


1997; 


Lery et al. 


2002) 



rials in this model are those deflected from infalling matter, instead of entrained or swept-up 
ambient gas by an underlying jet or wind. However, in at least one case, the G192.16 out- 
flow, there is strong evidence fo r the presence o f a wide-angle wind being the driving agent 
of the wide-angle ^^CO outflow (lOiu et al.ll2008l . and references therein). The 3.6 and 4.5 pm 



- 10 - 



emission coincident with the G240 ^^CO outflow lobes shown in Figure [2]i, is likely excited 
by an underlying fast wind. Therefore, we suggest that the G240 outflow is more likely to 
be driven/swept-up by an underlying wide-angle wind. The wind itself mi ght be accretion 



driven and somewh at analogous to X-wind or disk- wind in low-mass stars (IShu et al.l 12000 



Pudritz et al.l 120071 ). although it is unknown whether X-wind or disk-wind can be scaled up 



and applied to high-mass stars. 

Figure [S] shows the position- velocity (PV) diagram of the G240 ^^CO outflow cut roughly 
along the outflow major axis. The NW redshifted lobe shows a feature extending to ~ 5" 
from the driving source at a velocity of ~ 88 km s~^ and another feature extending to ~ 20" 
from the source at a velocity of ~ 80 km s^^. Connecting these two features following the 
gradient of the contour levels in Figure [H] one can get an approximately parabolic trajectory 
opening outward from the driving source. This kind of PV structure is very similar to that 
of the northern lobe of the VLA 05487 outflow, a good example of a wide-angle wind driven 
low-mass outflow, and ca n be reasonably reproduced by the hydrodynamic simulation of a 
modified X-wind model ( iLee et al.l 1200 ll ). According to the model the PV structure can 
be qualitatively understood in the context of a radially expanding shell with a Hubble-law 
velocity structure and a moderate inclination angle, where the far-side shell has a smaller 
projected extension and higher line-of-sight velocities, while the near-side shell has a larger 
projected extension and lower line-of-sight velocities. This parabolic PV structure supports 
the scenario that the G240 molecular outflow is driven by a wide-angle wind. It is difficult 
to identify a similar PV structure in the SE blueshifted lobe, probably because this lobe is 
highly asymmetric at most velocities (the western part of the shell can be only seen at 62 
and 64 km s~^ in Figure [3]) and contaminated by a secondary outflow to the west. Its PV 
structure is dominated by a high velocity feature ~ 14" away from the source, corresponding 
to a clump shown in Figure [2]d. 

While molecular outflows associated with both low- and high-mass (proto)stars are 
found to exhibit a "mass spectrum", m(f) ocf"'^, with the power-law index 7 around 2 for 
low velocity gas and often s teepening at veloc i ties higher than ~ 10 km s~^, its physical 
origin is not well established (IRicher et al.ll2000l : lArce et al.ll2007l ). Figure [9] shows the mass- 
velocity diagram of the G240 outflow, where both lobes can be fitted as a broken power-law, 
with 7 steepening from 4.3 ±0.1 to 7.0 ±0.3 at 13.5 km s~^ for the blueshifted lobe and from 
2.6 ±0.2 to 6.1 ±0.4 at 10.5 km s"^ for the redshifted lobe. It appears that the G240 outfiow 
has fairly large values of 7 at high velocities, comparabl e to that of some other ^^CO outfiows 



associated with both high- and low-rn ass (proto)stars (IShepherd et al.lll998l : IYu et al.ll2000 



Arce fc Goodmanll200ll : ISu et al.ll2004j ) . Such large values of 7 have been an unsolved issue to 



both jet bow-shock and wide-angle wind models since numerical simulations of these models 
predict 7 < 4; future simulations that run for fiows more than 10^ yr old might be able to 
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address whether it can be explained by outflow evolution (see a review by lArce et al.l 120071 ). 
Moreover, as shown by colored lines in Figure [9], the mass spectra of high velocity channels 
of both lobes can be further fitted with a broken power-law rather than a single power-law, 
and in those narrower velocity ranges 7 steepens from 5.9 ± 0.2 to 8.4 ± 0.7 at 19.5 km s~^ 
for the blueshifted lobe and from 5.3 ±0.3 to 9.6 ±0.8 at 20.5 km s~^ for the redshifted lobe. 
This may simply reflect a continuing decrease of mass entrainment efficiency with increasing 
outflow velocities in this outflow. 



4.3. Mass and energetics 



The total gas mass of the G240 outflow, estimated based on combined interferometer 
and single-dish data with sophisticated opacity correction, amounts to 98 M©, making it one 
of the most massive (proto) stellar outflows known. Given the unknown inclination angle of 
the outflow, the estimated outflow momentum of 620 M© km s~^ and energy of 4.5 x 10^^ 
erg can be treated as lower limits accordingly. Based on a velocity g radient in their CSO 
C^^O (2-1) observations and assuming the gas is in virial equilibrium, iKumar et al.l (120031 ) 
derived a gas mass of 625 Mq. Considering the cloud is likely infalling to feed the accretion 
onto the central (proto) stars and thus deviates from virial equilibrium, this estimate may 
provide a lower limit of the dense cloud ma ss. The gravitation al binding energy of the 
cloud, GM^i^^jR^ioud, where Rdoud = 0.63 pc JXumar et al.ll2003[ ). is then > 4.4 x 10^^ erg, 
which is comparable to the outflow energy (Table [1]). Although there are many uncertainties 
in calculations of both the gravitational binding energy and the outflow energy, the G240 
outflow seems to play a critical role in disrupting or unbinding its host cloud. The mass 
loss rate of 3.8 x 10"^ M© yr ^ is about two to three orders of magnitude higher than 
that of typical low-mass outflows. Assuming momentum conservation between the swept-up 
m olecular outflow and an underlying primary wind, i.e., Pout/idyn = M^Vw (also see Eq. 8 
of IShu et al.lll99ll ). where M^, are the wind mass loss rate and wind speed, respectively, 
and adopting a wind speed of 500 km s~^, the mass loss rate of the wind amounts to 
~ 0.5 X 10~^ Mq yr ~^. Furthermore, adopting a ra tio of | for the wind mass loss rate to 
mass accretion rate (ITomisaka Ill998l : IShu et al.ll2000l ). we get an accretion rate of 1.5 x 10~^ 
M© yr~^, high enough to overcome radiation pressure from the central object to form a 
high-mass star (jWolflre fc Cassinellilll987l : iJijina fc Adamslll996l : lYorke fc Sonnhalterll2002l ). 
Again considering the underestimate to Pout-, this estimate may represent a lower limit to 
the actual accretion rate. 
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5. Summciry 

We have presented a high-angular-resolution study of the molecular outflow in a 10"^'^ 
L0 star forming region, using the ^^CO and ^"^CO (2-1) observations from the SMA and 
CSO, as well as the SpitzerlRAC color-composite image. The SMA 1.3 mm dust continuum 
images are used to pinpoint the outflow powering source. 

A bipolar, wide-angle, quasi-parabolic molecular outflow is revealed in the ^^CO and 
^^CO emission at velocities more than 3 km s^^ apart from the cloud velocity. The outflow is 
most likely powered by a deeply embedded high-mass protostar prior to forming a significant 
UC Hii region. Prom its morphological similarity to some low-mass molecular outflows and 
parabohc PV structure in the redshifted lobe, the G240 outflow is hkely driven/swept-up 
by an underlying wide-angle wind. Por both lobes the outflow mass-velocity relation can 
be fltted with broken power-laws with the indices steeping at high velocities. The outflow 
has ~ 98 M0 gas mass and a dynamical timescale of 2.4 x 10^ yr, thus a mass loss rate of 
4.1 X 10^^ Mq yr^^. Without correction for the unknown inclination angle the momentum 
and energy of the outflow amount to 620 Mq km s~^ and 4.5 x 10"^^ erg, respectively. 

Our observations of the G240 outflow reveal a well-deflned wide-angle molecular outflow 
with its morphology and kinematics similar to low-mass molecular outflows, but with two 
to three orders of magnitude greater mass and energetics suggestive of a high accretion rate. 
These findings appear to be consistent with the scenario that high-mass stars of early-B 
to late-0 types can form in a similar way as low-mass star formation, even in a clustering 
environment. 

This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, Cahfornia Institute of Technology under a 
contract with the National Aeronautics and Space Administration (NASA). 
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07^44"^53.0^ 52.0^ 51.0^ 07"44"^52.20^ 52.00^ 51.80^ 

R.A. (J2000) 

Fig. 1. — : 1.3 mm continuum emission: (a) observed with the SMA compact array, contouring at 3x(-l, 1, 
2, 4, 7, 11, 16, 22, 29, 37, 46, 56, 67, 79, 92, 106) mJy bcam-\ where the r.m.s. noise level (Icr) is about 
1 mJy beam~^; (b) observed with the SMA extended array, contouring from 3.9 (3cr) to 54.6 mJy beam~^ 
at steps of 3.9 mJy beam~\ and the three mm continuum peaks are denoted as "MMl", "MM2", "MM3" 
from north to south. In the field all the negative contour levels are higher than — 3cr thus not shown. The 
corresponding synthesized beam is displayed in the lower left of each panel. 
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07"44"^53^ 51^ 
R.A. (J2000) 

Fig. 2. — : (a) low-velocity ^^CO 2-1 emission, integrated from 58 to 64 km s~^ for the blueshifted lobe (blue) 
and from 74 to 80 km s~^ for the redshifted lobe {red), with contours starting from 20% and continuing 
at steps of 10% of the peak emission, and two gray curves bisecting each other at the central mm peak 
outlines the shell of a bipolar wide-angle outflow, which is the focus of this paper; (b) high- velocity -'^^CO 
2-1 emission, integrated from 42 to 56 km s~^ for the blueshifted lobe (blue) and from 82 to 94 km s^^ 
for the redshifted lobe {red), with contours starting from 20% and continuing at steps of 10% of the peak 
emission, the gray contours show low- velocity ^^CO emission as shown in panel (a), and a straight gray line 
delineates the cut for a PV plot shown in Figure [51 (c) ^^CO 2-1 emission, integrated from 58 to 62 km s~^ 
for the blueshifted lobe {blue) and from 74 to 78 km s"""^ for the redshifted lobe {red), with contours starting 
from 15% and continuing at steps of 10% of the peak emission; (d) Spitzer IRAC color composite image with 
the 3.6, 4.5, and 8.0 /im emission coded in blue, greed, and red, respectively, the two curves are the same as 
shown in panel (a) to denote the molecular outflow shell, and the blue and red arrows mark the orientations 
of the blue- and redshfited lobes of the outflow, respectively. In (a) and (c) the corresponding synthesized 
beam is shown as an ellipse in the lower left. In each panel and hereafter, the central stars denote the three 
mm continuum peaks. 
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Fig. 3. — : ^^CO (2-1) channel maps from the SMA compact array observations. Sohd/dashed contours start 
from and continue at steps of 0.5/-0.5 Jy (+/-10cr) for channels of 40-58 km and of 78-94 km s^^, and 
1.0/-1.0 Jy for channels of 60-76 km s~^. Hereafter the central velocity of each channel is indicated in the 
upper left of each panel, and the synthesized beam is shown as an ellipse in the lower left of the first panel. 
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Fig. 4. — : Combined SMA and CSO ^^CO (2-1) channel maps of 60-76 km s ^, with sohd/dashed contours 
starting from and continuing at steps of 1.0/-1.0 Jy. 
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Fig. 5. — : The ^^CO (2-1) spectra at the CSO observational center, (R.A., decl.)j2ooo = 
(07''44'"52.1%-24'i07'"49''): the sohd black Hne shows the CSO spectrum; the soHd gray hne shows the 
spectrum constructed from the combined SMA and CSO data and convolved with the CSO beam; the dash- 
dotted gray line shows the spectrum obtained with the SMA compact array and convolved with the CSO 
beam. 
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Fig. 6. — : ^^CO (2-1) channel maps of 58-80 km s ^, from the SMA compact array observations. 
SoUd/dashed contours start from and continue at steps of 0.3/-0.3 Jy (+/-6tT). 
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Fig. 7. — : Combined SMA and CSO ^'^CO (2-1) channel maps of 64-70 km s ^, with soHd/dashed contours 
starting from and continuing at steps of 0.3/-0.3 Jy. 
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Fig. 8. — : Position- velocity diagram of the ^^CO outflow, cut along a position shown as a solid line in Figure 
[lb- The contour levels go from 1 to 22 Jy beam^^ in steps of 1.5 Jy beam^^, and the gray scale stretches in 
a square-root algorithm from 0.5 Jy beam~^ {white) to the peak emission at 27.8 Jy beam~^ (black). Plus 
signs outlines a parabolic trajectory (see the discussion in §4.21 for details). The vertical dashed line indicates 
Vcioud and the horizontal dashed line the position of the potential driving source of the outflow, i.e, the MM2 
peak. The diagram is constructed from the combined SMA and CSO data. 
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Fig. 9. — : Mass-velocity diagram of the "'^^CO outflow, calculated from the combined SMA and CSO data 
for the gas mass in each 2 km s~^ wide channel. Black filled squares denote the measurements from the 
blueshifted lobe, and gray filled circles the redshifted lobe. Solid black and gray lines indicate broken power- 
law fittings to the blue- and redshifted lobes, respectively; the blue and red lines represent further broken 
power-law fittings in the specific velocity ranges, rather than a single power-law fitting. 



Table 1. Derived core and outflow parameters 



source Menv^ 


M 2 






Mout 


Pout 


Eout 


tdyn 




(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo km s~i) 


(lO^^erg) 


(loV) 


(10-3 Moyr-^) 


G240.31+0.07 150 


40 


37 


61 


98 


620 


4.5 


2.4 


4.1 



^Mass of the dusty envelope as shown in Figure [IK 
^Combined mass of the three dusty cores as shown in Figure [T}d 
3 Outflow mass in the blueshifted lobe 
^Outflow mass in the redshifted lobe 



